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Abstract 

Within quantum mechanics model we study the problem of resummation of an 
asymptotic perturbation series for bound state parameters via optimization of the 
perturbative expansion. A possible application of the method to the positronium 
lifetime calculation is also briefly considered. 
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1 Introduction 



A discrepancy between the experimental result and theoretical predictions 0] for the 
orthopositronium width found quite a time ago still persists in spite of rigorous efforts 
to improve both theoretical computations and experimental data. By now an agreement 
on theoretical estimates for the width in the next-to-leading order is obtained and this 
agreed estimate is smaller than the experimental number by three standard deviations 
that causes some discomfort because it forces to allow quite big contribution of higher 
orders. The difference with experiment could be accounted for if future calculations of 
the order (a/V) 2 coefficient determine it to be of order 250 ± 40 [1J. The anomalously 
large next-to-leading correction can appear accidently in a sense that in higher orders 
corrections become small but it can also be a signal of bad convergence of the series in 
fine structure constant for the positronium bound energy (width). This does not seem 
impossible since in most of physically interesting models of quantum field theory the 
conventional perturbation theory forms the asymptotic series in coupling constant that 
can be used for calculation of the Green's functions only if the effective parameter of the 
expansion is small enough. However when the exact solution is absent it is difficult to 
determine if the asymptotic expansion is applicable to compute some physical quantity 
and the value a = 1/137 can be too large for the expansion in a being a good instrument 
in study of the orthopositronium width though it allows to compute another quantities, 
for example, the electron anomalous magnetic moment with high precision. Therefore it 
is instructive to try to go beyond the asymptotic expansion and improve the ordinary 
perturbation theory for orthopositronium width. 

Some methods have been suggested to improve a convergence of the conventional 
perturbation theory. The basic idea of the optimized S expansion is to introduce the 
artificial parameter 5 which interpolates between the theory we intend to solve with 
Hamiltonian H, and another theory, with Hamiltonian Hq(X) (A is a set of auxiliary 
parameters not present in the original theory), which is soluble and reflects the main 
properties of the theory we are interested in. One defines a new Hamiltonian depending 
on 5 
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H 5 = H (X) + 6(H-H {\)). 



(1) 



Then any desired quantity is evaluated as a perturbation series in 5, which is set equal 
to unity at the end of the calculations. The convergence of the series is achieved by 
optimization procedure i.e. by fixing the parameters A at every finite order of the 
expansion according to principle of minimal sensitivity (PMS) at the point where the 
result is least sensitive to their variation or principle of fastest apparent convergence 
(FAC) at the point where the next term in the series vanishes or somewhat else. Though 
above procedure is not rigorous it gives good numerical results in most of the cases. The 
method has been mostly advanced in studying the anharmonic oscillator [|[ where the 
convergence has been rigorously established jS|. 

Thus it seems instructive to optimize the perturbation theory in analysis of the positro- 
nium bound state. However the calculation of the correction to the orthopositronium 
width is very involved even in the case of ordinary a expansion. Therefore it is useful 
to consider the simplest model that, nevertheless, retains most relevant features of the 
real problem and, as we hope, can help to gain some intuition to cure the difficulties with 
positronium. 



The problem we will study is the ground state in spectrum of the stationary Schrodinger 
equation in three dimensions 



2 The model 




(2) 



where we imply 2m = % = 1 and 




r < 1, 



(3) 
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At r < 1 the potential (|j) consists of two parts: one is a constant and another depends 
linearly on r. The Schrodinger equation with the constant part of the potential only has 
bound states for any positive a with the ground state energy determined by the equation 



-E = yy + « cotU^- + a + E j . (4) 

Therefore if a is small enough we can search for the ground state energy of eq. (0) using 
perturbation theory and consider the constant part of the potential to be responsible for 
creation of the bound state while the linear term is a perturbation because it is suppressed 
by an extra power of a ( the parameter a should not be confused with the fine structure 
constant). 

On the other hand the exact solution of eq. (|2]) is known and leads to the equation for 
the ground state energy 

/ZW = 2/ 3/ Bi(^)Ai / fe)-Ai(e )Bi / te) \ 
^ Bi(£ )Ai(a) - Ai(£o)Bi(ei) J ' 

£o = -«~ 4/3 (j + a + E^j , £i = -a" 4 / 3 + a - a 2 + E^j (5) 

where Ai and Bi are Airy functions M . Using the asymptotic expansion of Airy functions 
at large negative (small a) we obtain an asymptotic series for the ground state energy 



, . ~. . a 2 ( /3 2 \ / 21 11 13\ o 
39 41 35 26\ , \ 

Substituting numerical values for the coefficients of the expansion (pi) 
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E(a) = -— (l - 1.7026a + 1.2602a 2 - 0.7864a 3 + . . .) (7) 

we find that the series merely reveals bad convergence near the point a ~ 1. In fact 
the series diverges for any positive a because the coefficients of the expansion (||) in high 
orders grow factorially. That reflects the presence of a singularity of the function E(a) 
at the origin of the complex a plane. The form of the singularity can be found directly 



4 



from eq. E(a) has a cut along negative semiaxis and a branching point at a = 0. 
For sufficiently small |a| it is an analytical function for — 7r < argct < tt therefore the 
series (|6]) is Borel recoverable i.e. we can extract complete information on the function 
E(a) from its asymptotic expansion. Presence of the singularity in the Green's function 
of eq. (0) reflects the fact that at a = the spectrum of eq. (Q) changes qualitatively and 
a discrete part of the spectrum appears. 

Such a singularity does not ultimately lead to divergence of the expansion of a bound 
state energy. For example, the ground state energy of the Schrodinger equation with 
the constant part of the potential (Q) is expanded in the convergent series for any finite 
positive a 

„ . , a 2 ( (\ 2 \ ( 5 11 3 \ „ 
E Q (a) = — - l- ---ra+ — - — - — a 2 - 



4 V2 tt 2 J ' V 16 6tt 2 tt 4 / 

/ 7 13 3 6 \ o \ 

-(32-8^-2^-^)" + ")- (8) 

So there is no implicit singularity in the series for Eq{q). However after inclusion of the 
perturbation such a singularity appears and the full series (|5]) becomes divergent. This, 
in a sense, simulates the positronium bound state where the non-relativistic Coulomb 
potential that is taken to build the leading order bound state Green's function does not 
lead to implicit singularity of bound state energy in the fine structure constant while the 
relativistic corrections result, as we suppose, in the divergent series. Thus we can use our 
quantum mechanics analog as a test model for further analysis of the positronium. 

3 Optimized expansion 

Our purpose now is to develop the optimized perturbation theory (OPT) for our toy 
model. Following the general idea we have to choose the form of the "unperturbed" 
Hamiltonian. The Hamiltonian 

Ho(a') = -A + U (a\ r) (9) 

where 
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, T + a') , r < 1, 
U (a', r) = { ^ 4 ' ' (10) 

0, r>l, 

is the potential of the spherical well with changeable depth seems to be the most appro- 
priate choice. Here the set of parameters A in eq. (|l|) is reduced to the single parameter 
a' characterizing the depth of the well. The ordinary perturbation theory corresponds to 
a 1 = a and 5 = 1. Making the expansion in 5 and setting 5 = 1 we obtain in n-th order 
the series for the ground state energy 

E n (a, a') = E^(a') + E^(a, «') + ••• + £ W («, a'), 



a' 2 / (\ 2 \ , / 5 11 3 



4 V V2 7T 2 / V 16 6tt 2 tt 4 



a' 2 



gW(a,aO = ^ [l-[;4)a' + -)+ (11) 
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(i + 4- (t-4 + ^) a + •••)- 



aV / 4 /3 2 12 



^WH(^+a*(i + i)) 2 + .... 

where the explicit expressions for are expanded in a' and a. In general if a' — a 
0(a 2 ) the effective parameter of the expansion is proportional to 



(0\U(a)-U o (a')\0) ^ 



^o(a') 

So E n (a, a') after expansion in a correctly reproduces the first n terms of eq. (^). Thus if 
we are interested in only the asymptotic expansion the choice of the start approximation 
does not play a crucial role. If, however, we intend to go beyond the asymptotic expansion 
we have to choose the "unperturbed" Hamiltonian to provide the best convergence of the 
expansion. The most transparent and conventional way to optimize the expansion ( |TT|) is 
to fix the parameter a' in n-th order according to PMS (FAC) at the value ot^ MS (<y^ AC ) 
so that 



^ =0, (13) 

oa a i= a PMs 
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or 



At n = 1 we have 



E^U 



PMS I i 

a, = a 1 



a FAC = a[l 



1 2 

2 + ^ ,a + 

1 2 

2 + ^'" 



Since a* — a = 0(a 2 ) where # stands for PMS or FAC we have E^ n \a, a*) = 
0{a n E^°\a)). Therefore expanding eq. (13, 14) in a we obtain at n — ► oo the 



series 



a 



# _y a # = a {\ + af a + af a 2 + 



moreover 



E n {a,a FAC )\ n ^^E Q {a FAC {a)). 

Now we show that one can choose a PMS = a FAC . Let us suppose that a PMS 
for i <n. The PMS condition to determine a p ^{ s reads 



On the other hand 



dE,, 



n+l 



da' 



a'-a PMS =0(a n + 3 ) 



O(a n+1 E (a)). 



E^\ a ^ aPMS=0{an+2) = O(a n+2 E (a)), 
E(a) — E n+ i\ al _ a PMs =0 ( a n+2- ) = O(a n+2 E (a)) 
because a PMS = a FAC for i < n. Therefore 



dR 



n+l 



da' 



d(E n+1 -E(a)) 



=0(a™+ 3 ) 



da' 



a'-a PMS =0{a n + 3 ) 



= 9(E n+ i\ a ^s - E(a)) + 0{a n + i Eoia)) = O (a^E (a)) (21) 

where taking a derivative of the asymptotic series is justified because E(a) is a power 
series and all non-analytical in a' terms in E n+ i must be cancelled in higher orders and, 
therefore, can be omitted before taking a derivative. So condition (|19"D is satisfied for 
arbitrary a™[ s and we can put a^[ s = a^p. Since af MS = af AC for i = 1 (see 
eqs. (|15[ fL6|)) one can choose af MS = af AC for any i whence a PMS = ot FAC '. 

Thus in the limit n —>■ oo all corrections in the optimized expansion ( |TT| ) vanish both 
for PMS and FAC optimization prescription and we obtain 

E n (a, a*) -> E (a # (a)) = £7(a). (22) 

Since the function E Q (a#) can be expanded in the convergent series in a # (eq. (§)) 
the singularity of the function i£(a) at a = is absorbed by the function or (a) i.e. the 
series fll6|) is a divergent asymptotic expansion. However we can search for the values 
a* which satisfy eqs. ( |I3| , P^ ) at every order numerically rather than as a series in a. In 
Fig. 1 and Fig. 2 we plot the functions ^(1.0, a!) and E^(l,a')/E^(l,a') to show a 
typical pictures that are used to determine a # . Thus for a = 1.0 we obtain af M5 = 0.310 
and a£ AC = 0.304. 

Thus we have three kinds of perturbative expansion: the asymptotic series, the stan- 
dard perturbation theory and OPT. The asymptotic series seems to be the most primitive 
tool and we expect to obtain the best results using OPT. This assumption is completely 
confirmed by numerical analysis given in the next section. 



4 Numerical evaluation 

The results of numerical analysis are given in Tabs. 1 — 30. In Tab. 1 we establish the 

exact value E(a) along with the results of the asymptotic expansion E n (a) up to the n-th 

order (n = 0, 1, 2, 3). In Tab. 2 we compare the exact value E(a), the result of the 
1 For the exact value and approximations of the bound energy in Tabs. 1 — 3 the factor (— 10~ 2 ) is 
implied 
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optimized expansion Ei(a,a PM ) and the results of the odinary perturbation theory up 
to the n-th order (n = 0, 1) i.e. not expanded values of E (a) and Ei(a, a). All the [i,j] 
Pade approximants E^'^(a) with i + j < 3 are in Tab. 3. 

We give numerical estimates for three different values of the parameter a that represent 
typical cases 

1. a = 0.1: the asymptotic expansion is applicable and justified; 

2. a = 0.5: the asymptotic expansion is applicable but one has to deal with the high 
orders corrections to achieve a satisfactory accuracy. An improvement of the perturbation 
theory is desirable; 

3. a = 1.0: the asymptotic expansion in principle can provides only ~ 2% accuracy after 
summation of ~ 13 terms. Then the terms start to grow. The perturbation theory must 
be reformulated. 

As we can see a naive attempt to improve the convergence of the perturbative expan- 
sion simply keeping the exact (not expanded) value of E n (a, a) that sums up some next- 
to-leading corrections gives a good result but is essentially insufficient if a is large enough. 
This shows that optimization is important for convergence. Though we cannot directly 
demonstrate that without the optimization the series (|TT|) becomes divergent as it takes 
place in the case of anharmonic oscillator || the numerical analysis clearly shows an ad- 
vantage of optimized expansion. Indeed in all these cases the best convergence is achieved 
within OPT. Even for a = 1.0 taking only the first order correction we reach ~ 2% ac- 
curacy. Using the PMS prescription we can also correctly estimate an error of the result. 
For example for a = 1.0 the naive estimate is (£ (1) (1.0, a PMS )/E^(1.0, a PMS )) 2 ~ 0.02 
that coincides with the real uncertainty (see Tab. 1). 

We should note that in all the cases the terms of a expansion taken into account in 
our numerical analysis are far from the critical point where the series begins to diverge. 
The bad convergence of the series reveals only in the fact that they decrease quite slow. 
On the other hand even the third order correction is hardly available within ordinary 
perturbation theory. So the accuracy of the perturbation theory is restricted rather by 
technical reasons than by asymptotic character of the series. Thus the optimization is not 
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only a resummation prescription that is useful in high orders of asymptotic expansion but 
also gives an opportunity to improve accuracy of perturbation theory in low orders. This 
is an important benefit especially for non-trivial systems where high orders calculations 
are impossible. 

A remark about Pade approach is in order. As we can see some Pade approximants 
are closer to the exact result than the plain asymptotic expansion. However it is not 
possible to make a choice between various approximants until the exact result or the 
general structure of the series are known. Moreover in high orders where the asymptotic 
character of the series reveals the Pade theory becomes useless. The matter is that Pade 
approximants because of their specific structure can correctly reproduce only pole-like 
singularity while the function E(a) has a branching point at a = 0. 

5 The positronium bound state 

The theoretical predictions for the orthopositronium width reads [0 



where the coefficient B has not yet been computed but it is expected to be about 250 to 
bring theory and experiment into agreement. The coefficients of eq. (|23|) tend to grow 
rapidly. In spite of the value a = 1/137 seems to be small enough the contribution of the 
second order term is about 0.2% of the leading one. This resembles the behavior of the 
asymptotic series in our toy model with 0.1 < a < 0.5 where the asymptotic expansion 
can, in principle, provide a sufficient accuracy but it requires to compute high order terms 
that is not possible by pure technical reasons. So we hope that after optimization of the 
expansion an agreement between experimental and theoretical values will be reached 
already in the first order of OPT. 

The main problem of OPT is to find an appropriate form of "unperturbed" action 
which, on the one hand reflects the main properties of the exact theory and, on the other 





(23) 
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hand does not lead to extremely cumbersome calculations. The most direct way for the 
positronium bound state is as follows. By now there are two general methods to develop 
the systematic perturbation theory for the positronium: the non-relativistic 1/c expansion 
H and the method based on Bethe-Salpeter equation ||. The non-relativistic Coulomb 
solution and Barbieri-Remiddi solution [JXD | which are used to built the leading order 
approximations for the positronium bound state Green's function within these approaches 
involve the physical fine structure constant as a parameter. In analogy with our toy 
model one can replace it with a changeable parameter a'. Then in every order of the 
new perturbation theory one has to keep the exact dependence on a' and fix it using 
some optimization prescription. We should emphasize that the optimized value has no 
direct relation to the physical fine structure constant. This is an auxiliary parameter and 
constructing the OPT, for example, for orthopositronium bound energy we would obtain a 
different value of this parameter. Unlike the case of the asymptotic a expansion two above 
start approximations can lead to the different results if one is interested in the optimized 
expansion but it is not obvious which is preferable. Though the detail analysis of the 
positronium is a non-trivial technical problem that is a subject of a separate publication. 
To conclude we should note that an attempt to improve the perturbation theory in 



orthopositronium width analysis using Pade approximants has been made in ref. [11 



However it looks artificially because an information about structure of the asymptotic 
expansion ( p^ ) is absent. On the other hand OPT seems to be the most appropriate tool 
to deal with such a problem because it speeds up the convergence of the perturbation series 
choosing the most natural start approximation for a specific model not by an extraneous 
mathematical trick. This is an automatic summation device that need not an input 
information on the form of the singularity of the bound state Green's function in the 
coupling constant but reproduce it via optimization. This feature can be merely observed 
in our toy model where the form of the singularity of the ground state energy at a = in 
the auxiliary and original theories are essentially different but OPT reproduces the correct 
singular a dependence that is absorbed by the optimized value of auxiliary parameter a'. 
This general property of optimized perturbation theory allows to cope even with Borel 
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non-summable series || while a class of the problems where the Pade theory can be 
successfully applied is quite restricted. 
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Figure Captions 

Fig. 1. The curve line depicts the function Ex(1.0, a'). The straight line corresponds to 
the exact value £7(1.0). 

Fig. 2. The function E^(1.0,a')/E^(1.0,a'). 
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Tables 



Table 1. The exact bound energy E(a) and the results of the asymptotic expansion up 
to the n-th order E n (a) (n = 0, 1, 2, 3). 

Table 1. 



a 


E(a) 


E (a) 


Ma) 


E 2 (a) 


E 3 (a) 


0.1 


0.2101 


0.2500 


0.2074 


0.2105 


0.2104 


0.5 


2.426 


6.250 


0.9292 


2.898 


2.284 


1.0 


2.144 


25.00 


-17.57 


13.94 


-5.721 



Table 2. The exact bound energy E(a), the result of the optimized expansion up to the 
first order E\ and the results of the odinary perturbation theory up to the n-th 

order E n (a, a) (n = 0, 1). 



Table 2. 



a 


E(a) 


E^a™ 3 ) 


E (a) 


Ei(a, a) 


0.1 


0.2101 


0.2104 


0.2430 


0.2092 


0.5 


2.426 


2.424 


5.531 


1.915 


1.0 


2.144 


2.116 


20.13 


-4.200 



Table 3. The exact bound energy E(a) and various Pade approximants E^\a) with 
i+j <3. 

Table 3. 



a 


E(a) 


E^(a) 


E^ 2 \a) 


E^ l \a) 


$ (a) 


E^(a) 


0.1 


0.2101 


0.2104 


0.1874 


0.2104 


0.2104 


0.2104 


0.5 


2.426 


2.366 


2.340 


2.430 


2.424 


2.562 


1.0 


2.144 


0.5388 


5.759 


1.833 


1.600 


4.331 
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Fig. 2 



